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If vector type perturbations are present in the primordial plasma before recombination, the gen¬ 
eration of magnetic fields is known to be inevitable through the Harrison mechanism. In the context 
of the standard cosmological perturbation theory, nonlinear couplings of first-order scalar perturba¬ 
tions create second-order vector perturbations, which generate magnetic fields. Here we reinvestigate 
the generation of magnetic fields at second-order in cosmological perturbations on the basis of our 
previous study, and extend it by newly taking into account the time evolution of purely second-order 
vector perturbations with a newly developed second-order Boltzmann code. We confirm that the 
amplitude of magnetic fields from the product-terms of the first-order scalar modes is consistent 
with the result in our previous study. However, we find, both numerically and analytically, that the 
magnetic fields from the purely second-order vector perturbations partially cancel out the magnetic 
fields from one of the product-terms of the first-order scalar modes, in the tight coupling regime 
in the radiation dominated era. Therefore, the amplitude of the magnetic fields on small scales, 

^ ^ 10 /iMpc“^, is smaller than the previous estimates. The amplitude of the generated magnetic 
fields at cosmological recombination is about Brec = 5.0 x Gauss on fc = 5.0 x 10~^ ZiMpc”^. 

Finally, we discuss the reason for the discrepancies that exist in estimates of the amplitude of 
magnetic fields among other authors. 


I. INTRODUCTION 

The presence of magnetic fields on large scales is established by current observations [H-Q. Such cosmological 
magnetic fields coevolve with the Universe, e.g., astrophysical objects, cosmic microwave background radiation (CMB), 
large scale structure, and inflation. Recent observations indicate that cosmological magnetic fields have the strength 
about micro-Gauss on Mpc scales (see, e.g., Refs. [la&i , and references therein). Moreover, the pair-echo method 
asi determines the lower bound of strength in the intergalactic magnetic fields a.s B > 0(10 Gauss. The 
remarkable progress of observations indicates that cosmological magnetic fields appear everywhere even in the cosmic 
voids. Furthermore, the future experiments of the radio telescope such as Square Kilometer Array can survey much 
deeper and wider regions, and give us rich information about coevolution between cosmological magnetic fields and 
baryonic matters. However, very interestingly, the origin of cosmological magnetic fields is not entirely revealed. 

The key process related to evolution of magnetic fields is the dynamo mechanism [3,[I34I3) which is the amplification 
mechanism of magnetic fields in the nonlinear magnetohydrodynamics. In the stars, galaxies, and galaxy clusters, 
their nonlinear evolution can amplify seed magnetic fields, and the strength of seed fields is about 10“^° ~ 10“^° 
Gauss [i3- Generally, the dynamo mechanism cannot generate magnetic fields from the absence of seed fields but 
amplify the seed fields, which should be set before the dynamo mechanism works. Therefore, when we believe that 
the origin of cosmological magnetic fields is as a result of amplification of seed fields by the dynamo mechanism, seed 
fields must be created in the early stage of the universe, namely before cosmological recombination. 

One of the candidates to create seed fields is the quantum fluctuations of the electromagnetic fields in the inflation 
era. During the inflation era, the scale of fluctuations is extended beyond the Hubble horizon due to the nearly 
exponential expansion. At first glance, it is possible to rely on this scenario to generate large-scale magnetic fields. 
However, this scenario does not work since the standard Maxwell theory has the symmetry under the conformal 
transformation. Under this symmetry, the vector fields, such as magnetic fields, undergo decaying only and become 
negligible. In other words, when the conformal invariance is broken, inflationary magnetogenesis possibly works. In 
many of the previous studies (Tsl - f^ , the authors have introduced interaction between the electromagnetic fields and 
the dilaton-like scalar field to break conformal invariance. However, even in this case, there are other problems in 
the inflationary magnetogenesis, i.e., strong-coupling and backreaction problems [H, These problems make the 


Electronic address: saga.shohei@nagoya-u.jp 



2 


situation worse. Therefore we can conclude that it is difficult to generate seed fields during the inflation era alone. 
In the more recent studies [25l - l^ , a new interaction between electromagnetic fields and axiom like pseudoscalar field 
is added in the context of dilaton-like magnetogenesis. This interaction ends up with generation of helical magnetic 
fields on small scales. The inverse cascade can transfer helical magnetic fields from smaller scales to much larger scales 
owing to conservation of the magnetic helicity. We therefore expect the existence of magnetic fields on all scales in 
later epochs although the detailed numerical estimation is needed [H, [ 2 ^. 

The cosmological phase transition is another possibility to generate seed fields (e.g., Refs. [3, 0). In general, 
phase transitions release the free energy and electric charges. The released free energy is converted into the electric 
currents. If these electric currents have rotational components, cosmological seed fields are induced at epochs of phase 
transitions. However, the coherent length of seed fields generated in the cosmological phase transitions cannot exceed 
the Hubble horizon scale at that time due to the causality. Therefore phase transitions alone are not able to explain 
observed large-scale magnetic fields [s^. 

Another category of generation mechanism is originated from astrophysical phenomena. For example, the Biermann 
battery is one of the candidates of generation mechanism after recombination. The gravitational force can be described 
as a gradient of the scalar potential and hence cannot generate vorticity. However, the Biermann battery, which 
is nonadiabatic phenomena such as shocks, can generate vorticity and subsequently, seed helds are induced with 
the amplitude about 10“^^ Gauss in protogalaxies [3l|, 10“^^ ^ 10“^^ Gauss in supernova remnants [s^, and 
~ 10“^^ Gauss in galaxies [S^. The Weibel instabilit y, w hich is microscopic instability in the plasma, can amplify 
tiny seed fields at the epoch of structure formation . When the velocity distribution of plasma particles has an 

anisotropy in the phase space, the isotropized process of the velocity distribution releases the energy, and subsequently, 
the energy is converted into magnetic fields. Accordingly, magnetic helds amplihed by the Weibel instability have a 
quite large amplitude of about 10“^ Gauss [S^, [1^. However, the Biermann battery and Weibel instability can only 
work with the existence of baryonic matters or astrophysical objects. Therefore it is difficult to explain the origin of 
intergalactic magnetic helds or magnetic helds in the voids. 

Yet another interesting mechanism of generating magnetic helds is the Harrison mechanism in which magnetic 
helds are generated via vorticity of the primordial plasma. In Ref. [^ . the authors have formulated the Harrison 
mechanism based on the cosmological perturbation theory in the primordial plasma which is a multicomponent 
system composed of photons, electrons, protons, dark matters, and neutrinos. In this system, photons and electrons 
or protons interact with each other through the Compton scattering. However, photons push electrons more frequently 
than protons because of the difference of scattering rates. Accordingly, the charge separation takes place. If there exist 
rotation-type electric fields, magnetic fields are generated. However, in the linear perturbation theory, the Harrison 
mechanism does not work because there is no growing mode solution for the vector-mode perturbations which induce 
rotation-type electric fields. In other words, the models including the active vector mode supplied by external sources, 
i.e., free-streaming neutrinos [38l-l40l|. cosmic defects [il|,|43, and modified gravity with vector fields [i^, can generate 
magnetic fields via the Harrison mechanism. 

Moreover, it turns out that even standard cosmological perturbations can generate magnetic fields if we take into 
account contributions from higher-order perturbations. In fact, it is known that the second-order perturbation theory 
has not only the scalar mode but also the vector and tensor modes through the product of the first-order scalar 
perturbations. Recently, the second-order cosmological perturbation theory is well established in the context of the 
GMB formalisms |43 - l^ . For example, the B-mode polarization is calculated based on the second-order perturbation 
theory while there is no B-mode polarization in linear scalar perturbations. 

Recently, generation of magnetic fields via the second-order perturbation has been studied in detail In 

these studies, the tight-coupling approximation is employed to estimate the amplitude of magnetic fields analytically. 
Each study has shown that the amplitude of generated magnetic fields is about 10“^° ^ 10“^^ Gauss at recombination 
on Mpc scales. However, it is difficult to know the detail of the magnetic power spectrum since the tight-coupling 
approximation breaks down inside the horizon scale at recombination. By solving perturbation equations up to the 
second order without employing the tight-coupling approximation, it is possible to analyze the power spectrum of 
magnetic fields. In Ref. [^, the authors have evaluated the spectrum generated by the vorticity of charged particles, 
which are induced by the nonlinear coupling between the first-order density perturbations. And they have found 
that resultant comoving magnetic fields have the amplitude of about 10“^® Gauss at recombination on Mpc scales. 
Subsequently, in Refs. |62l. l63l|. the authors have studied the Harrison mechanism including the anisotropic stress of 
photons. They have found that the amplitude of magnetic helds has 10“^° Gauss at recombination on Mpc scales. 
However, they ignore the purely second-order velocity difference between charged particles and photons in their 
analysis. In Ref. [b^, the authors include the purely second-order effects for the hrst time and analyze the spectrum 
of magnetic helds on superhorizon scales. In these studies, however, there are some discrepancies which have to be 
clarihed. 

In this paper we numerically solve the vector mode of cosmological Einstein-Boltzmann equations at the second 
order including all the effects relevant to the generation of magnetic helds, with a newly developed numerical code. In 
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addition we present analytic interpretations of the shapes and time evolutions of the power spectrum of magnetic helds 
on sub- and superhorizon scales, and make it clear what has caused the discrepancies among the previous studies. 
The paper is organized as follows. In Sec. II, we review the generation mechanism of magnetic fields. We decompose 
the evolution equation of magnetic fields up to the second order in terms of scalar and vector perturbations. It 
is shown that the magnetic fields are generated by the vector mode only. In Sec. Ill, we describe the perturbed 
Einstein-Boltzmann system up to the second order to compute the purely second-order perturbations, and show the 
tight-coupling solutions of the system. Time evolutions of magnetic fields and spectra are numerically evaluated and 
the analytic expressions are given in Sec. IV. We devote Sec. V to our discussions and conclusions. 

Throughout this paper, we use the units in which c = h = 1 and the metric signature as (—,-I-,-|-,-I-). We obey 
the rule that the subscripts and superscripts of the Greek characters and alphabets run from 0 to 3 and from I to 3, 
respectively. 


II. GENERATION OF MAGNETIC FIELDS 

In this section, we review basic equations for the generation of magnetic fields [^1^. i.e., perturbation equations of 
photon, proton, and electron fluids. While protons and electrons are conventionally treated as a single fluid, however, 
it is necessary to deal with proton and electron fluids separately in order to discuss the generation of magnetic fields. 
Let us begin with the Euler equations. Those are given by 

mpnu^Upi-f^ - enu'^Fip, = -I- (1) 

nienu^Uei-fi + enu'^Fip, = C'®f -|- C®/ , (2) 

where mp(e) is the proton (electron) mass, Wp(e) is the bulk velocity of protons (electrons), Fp^i is the usual Maxwell 
tensor. The thermal pressures of proton and electron fluids are neglected. The right-hand side of Eqs. (ED and (ED 
represent the collision terms. The first terms in Eqs . ED and ([2]) are collision terms for the Coulomb scattering 
between protons and electrons, which are given by [56l| 

Cf® = -(7®^ = -{ui- Uei)e^n^ri, , (3) 


where 


riv = 


1 1/2 
•Ke me! 


In A 


9.4 X 10 ^®sec 




( 4 ) 


is the resistivity of the plasma and InA ~ 0{1) is the Coulomb logarithm, which is the almost constant parameter. 
As is well known, this term acts as the diffusion term in the evolution equation of magnetic field. The importance of 
the diffusion effect can be estimated by the diffusion scale, 


Adiff = \/rvr ^ 100 



( 5 ) 


above which magnetic field cannot diffuse in the time-scale r. Here FIq = lOO/i km/s/Mpc is the present Hubble 
parameter with h being the normalized Hubble parameter. Thus, at cosmological scales considered in this paper, this 
term can be safely neglected. 

The other terms expressed by are the collision terms for Compton scattering of protons (electrons) with 

photons. Since photons scatter off electrons preferentially compared with protons by a factor of (me/mp)^, we can 
safely drop the term from the Euler equation of protons. This difference in collision terms between protons and 
electrons ensures that small difference in velocity between protons and electrons, that is, electric current, is indeed 
generated once the Compton scattering becomes effective. In the next subsection, we derive the explicit form of the 
Compton scattering term. 


A. Gompton collision term 

Let us now evaluate the Compton scattering term. In the limit of completely elastic collisions between photons 
and electrons, this term vanishes. Typically, in the regime of interest in this paper, very little energy is transferred 
between electrons and photons in Compton scatterings. Therefore, it is a good approximation to expand the collision 
term systematically in powers of the energy transfer. 
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Let us demonstrate this specifically. We consider the collision process 

7(p^) + e"(g^)7(p'^) + (6) 

where the quantities in the parentheses denote the particle momenta. To calculate this process, we evaluate the 
collision term in the Boltzmann equation of photons: 


C[f] 


d^p' 


d?q 


d^q' 


E{p) J {2^f2E{p') {27r)32Ee{q) {2^f2E,{q') 


\Mr 


x{2-Kf5'^{q>^ +p>^ - q'^ -p''") ffe(g')/(p')(l + fiP)) 


5e(g)/(p)(l + /(p')) 


( 7 ) 


where f{p) and ge(g) a-re the distribution functions of photons and electrons, Ee(q) = + is the energy of 

an electron, and the delta functions enforce the energy and momentum conservations. And |A1|^ is the scattering 
amplitude for Compton scattering. We have dropped the Pauli blocking factor (1 — pe). The Pauli blocking factor 
can be always omitted safely in the epoch of interest, because pe is very small after electron-positron annihilations. 

Integrating over q', we obtain 


^ -pj 

xS 


d^p' 


d^q 


27r 


\M\^ 


{2Tr)^2p' (27r)32£'e(g) 2Ae(|q + p - p‘ 
p-p' + Ee{q) - Ee{\q +p- p'\) \ \peiq +P- p')/(p')(l + fiP)) - 5e(g)/(p)(l + /(p')) 


( 8 ) 


In the regime of our interest, energy transfer through the Compton scattering is small and can be ignored in the first 
order density perturbations. As we already discussed earlier, however, it is essential to take the second-order couplings 
in the Compton scattering term into consideration for the generation of magnetic fields. Therefore, we expand the 
collision term up to the first order in powers of the energy transfer^, and keep terms up to the second order in density 
perturbations. 

The expansion parameter is the energy transfer. 


Ee{q) - Eei\q + p - p'\) 


jp' -p) q _ (p - p') 

nrie 2me 


(9) 


over the temperature of the universe. Employing p T, we can estimate the order of this expansion parameter as 

~ 0{^)^ which is small when electrons are nonrelativistic. Note that, in the cosmological Thomson regime, 

2 

electrons in the thermal bath of photons are nonrelativistic, p ~ 2 ^’ energy of photons is much smaller than 

the rest mass of a electron, p <C mg. Thus, it also holds that q ^ ^/2rn^ p, and the second term in Eq. (|9]) is 
usually smaller than the first one. 

Now let us divide the collision integral into four parts, i.e., the denominators of the Lorentz invariant volume, 
the scattering amplitude, the delta function, and the distribution functions, and expand them due to the expansion 
parameter defined above. First of all, the denominator in the Lorentz invariant volume can be expanded to 


where 


1 

Ee{q)Ee{\q + p - p'\) 




2mt 


\q + p 


1 

2 

mi 


1 — 







(p -p')q 

9 

mi 



{P-P'f 

2mi 


( 10 ) 


( 11 ) 


Second, we consider the scattering amplitude. Fortunately, it has been known that the leading term (zeroth order 
term), obtained by multiplying together the first term in the delta function and the zeroth-order distribution functions, 
is zero. It means that we only have to keep up to the first order terms when we expand the scattering amplitude and 


^ However, we shall keep up to the second-order terms for the purpose of deriving the Einstein-Boltzmann system in section [m] 
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the energies, in order to keep the collision term up to the second order [66l| . The scattering amplitude for Compton 
scattering in the rest frame of the electron is given by, 


= &TTm\(jT 


cos p = p ■ p , 


P , P ■ 1 a 
-h — sm a 

P p' 


( 12 ) 


where p and p' are the energies of incident and scattered photons, p and p' are the unit vectors of p and p', respectively, 
denoting the directions of the photons in this frame. The Lorentz transformation with electron’s velocity (q/mg) gives 
the following relations. 


P ^ x/1 - {q/mgY 

P l-p-q/(pTOe) ’ 


P^P^ = P^Pti ■ 

Using these relations, we evaluate the scattering amplitude in the CMB frame as [b? 

= G-TTmeCrr Mo+M(^^) , 

where 


= 1 + cos^/3 , A4( 1 ) = —2cosd(l — cos/3) 

V me ^ 

Here ft and n' are the unit vectors of p and p', respectively. 

Third, we expand the delta function to 


— • (ri + n') 

rUe 


Hp-P'+ Ee{q) - Ee{q')] « 6{p-p') + 


dS [p-p' + Ee{q) - Eg{q')] 


dp 


{Eg{q) - Eg{q')) 


q=q' 


1 d'^S [p-p' + Egjq) - Eg{q')] 

2 dp^ 


{Eg{q) - Eg{q')f 


d5(p-p'). 


dd{p — p') 


= S{p - p') + " D(_3_) + """U . " D(^^ + -V\ 


d'^5{p — p') 


where 


(—) 
' TTLo ' 


dp' 

[p -p')q 


V, 


^ dp' ^ 2 

(P - P')^ 


( —) 
'• me. ' 


2me 


Finally, the distribution of the electron can be expanded to 


9e{q+p- p') ~ ffe(g) + ^ • (p - p') + ^ (P* - P'^^ QqiQqj ~ P^)' 

We assume that the electrons are kept in thermal equilibrium and in the Boltzmann distribution: 


9e{q) = ' 


f—) 

\meTeJ 


3/2 


exp 


(g - UleVg) 


21 


2meTe 


(13) 

(14) 

(15) 

(16) 


(17) 


(18) 


(19) 


( 20 ) 


where Vg is the bulk velocity of electrons. The derivatives of the distribution function with respect to the momentum 
are given as 


dqe _ _ g» - rrigVgi 

dq' rUgTg 

d'^Qe _ dge q' - rngvl 5'^ 

dq'dqi dqi rrigTg rrigTg 


( 21 ) 


( 22 ) 
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By substituting the above equations, Eq. dl) is written as 


ge{q + p-p') ~ ge{q) 


1 - 

' ' 2 ^ '' ^ 


where 


Therefore, we have 


q- rrieVe , 1 (p - p' f 

^{^) = - 7^ - (P-P), -F(^) = -- - - 

rriele ® 2 mgie 


5e(g + P - p')/(p')(l + /(p)) - 5e(g)/(p)(l + /(p')) 

= 5e(9) [/(p') - /(p)] - /(p')3e(g)^(^) - /(p')5e(g) 

1 


2 me '' 


+5e(g)/(p)/(p') 




( 23 ) 


(24) 


(25) 


Combining altogether, we obtain the collision term expanded with respect to the energy transfer as (note that this 
expansion is not with respect to the density perturbations) 

C[f] = ^TT^— [ , f [(0th order) + (1st order) + (2nd order)] , (26) 

2 p J (27r)V J {2ttY 

where 


0th order term; 


1st order terms: 


2nd order terms: 


Mo5{p - p')ge{q) [/(p') - /(p)], 


Mogeiq) 


-Sip - p')/(p')-^(^) + ^ ^ [fip') - fiP)] 'D(^) 


+M^^^geiq)Sip - p) [fip') - fip )], 


Mogeiq) 


-Sip-p')fip') ifip') - fip)] 


, d6ip-p')^ dSip-p')^ 


+7W(^)ge(g) 
+Mogeiq)fip)fip') 


-Sip - p')/(p')-7^(^) + 


dSip — p') 


-geiq)fip)fip')M(^^)J^(^)Sip - p). 


(27) 


(28) 


(29) 


From now on, we omit the second-order terms. These terms are not only much smaller than the first-order terms 
but also may not contribute to the Euler equation at all (see. Ref. El)-' Evaluating the first moment of the above 
collision term, we obtain the Compton scattering term in the Euler equation (1371) as 


cr = - 


[ d^P 

J (2^)3 

4aneaT 


P^C[f] 


-Py 


3 

('^ez ’^7^ 


(30) 


where the product of the velocity of electrons and anisotropic stress of photons in Eq. (1501) is the anisotropic part 
of “radiation drag” in the context of the radiation hydrodynamics. The radiation drag is originated by the electron 
motion in anisotropic radiation fields with absorptions and emissions. The velocity of electrons Vei obeys the Euler 
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equation given in the next section, and the velocity and anisotropic stress of photons, v^i and , obey the Boltzmann 
equation. We will show the explicit equations for electrons and photons in Sec. IIII Al And we will show the anisotropic 
stress of photons can be written in terms of the brightness function, which is also defined in Sec. IIII Al 
Here moments of the distribution functions are given by 


Pfj{p) = Pi, 


^ Pi flip) — nPl^li, 


dfp 

w 

dfp 
(271-) 

(Pq 

Pq 

P PiPj flip) ~ Pi^iij 


geiq) = Ue, 
^qiPeig) = PeVe 


(271-) 


“b 2 Pi^ij j 


(31) 

(32) 

(33) 

(34) 

(35) 


where pj and Pe(= rneUe) are energy densities of photons and electrons, vl^ and u® are their bulk three velocities 
defined by v'^ = u'' jp, and Hlj? is anisotropic stress of photons. Hereafter we use the relative velocity between photons 
and electrons as — Vei- It should be noted that the collision term (1501) was obtained nonperturbatively 

with respect to density perturbations 


B. Evolution equations of magnetic fields 


Now we obtain the Euler equations for protons and electrons as 


,, „ „ AaTPian 

TTienu^Uei-p + enu^Pif^ = --- 




(36) 

(37) 


where rup is the proton mass. Here, we ignore the pressures of proton and electron fluids. In addition, the Coulomb 
collision term is neglected as explained below Eq. ®. Note that the collision term was not evaluated in a manifestly 
covariant way. Here the left-hand side in Eqs. (|M)) and (1^ should be evaluated in a conformal coordinate system. 
We also assumed the local charge neutrality: n = Ue ^ Up. In the case without electromagnetic fields (Fi^ = 0), the 
sum of the equations (13611 and (13711 gives the Euler equation for the baryons in the standard perturbation theory. On 
the other hand, subtracting Eq. (15511 multiplied by me from Eq. (15711 multiplied by rUp, we obtain 


np ( — 
n 








rrip - rrie Ji 

mp H- rUe en 


+ en(mp -|- me)PF^^ - imp - me)j^Fi 


- Ifl 


Arripp^anaT 




(38) 


where P and are the center-of-mass 4-velocity of the proton and electron fluids and the net electric current, 
respectively, defined as 


mpP + meU^ 

P = --, 

mp + me 
f ^ en{Pp - <). 


(39) 

(40) 


Employing the Maxwell equations F ^''we see that the quantities in the square bracket in the left-hand side 
of Eq. (1551l is suppressed at the recombination epoch, compared to the second term, by a factor 


Pul 


3x10-40 


(41) 


where c is the speed of light, A is a characteristic length of the system and Wp = \fFKnP^Jme is the plasma frequency. 


























The third term in the left-hand side of Eq. (l38l) . i.e., {nip — me)j^Fip_, is the Hall term which can also be neglected 
because the Coulomb coupling between protons and electrons is so tight that |m*| » [Up — itg|. Then we obtain a 
generalized Ohm’s law: 




Aarp-yCi 

3e 


6V-y\yi , UptH, 




= a. 


(42) 


Now we derive the evolution equation for the magnetic field, which can be obtained from the Bianchi identities 
= 0, as 


0 = 

= - .«* U., + 




(43) 


where is the Levi-Civita tensor and H® = = e^^^Fjk/2 is the magnetic field in the comoving frame [^. We 

will now expand the photon energy density, fluid velocities, and photon anisotropic stress with respect to the density 
perturbation as 


p^{t,Xi) = -, u°{t,Xi) = a{t) ^-I-u(^^°(t,a;i)-I-, 

u'^{t,Xi) = u^'^'>’-{t,Xi) + ^u^'^'>’'{t,Xi) H-, Vi{t,Xi) = vl^\t,Xi) + i?;f^(t,Xj) H- 

ni;^'(t,x0 = nW®^'(t,x,) + -- -, (44) 

where the superscripts (0), (1), and (2) denote the order of expansion and t is the cosmic time. Remembering that H® 
is a second-order quantity, we see that all terms involving H® in Eq. (l4^ . other than the first term, can be neglected. 
Thus we obtain 


dg® 

dt 


««* c„ + 


/I (0) 

3e 




i(' 


el 


(45) 


where we used the density contrast of photons, = p^'l/p^'^- Eurther, we employed the fact that there is no 
vorticity in the linear order: = 0. It should be noted that the velocity of electron fluid can be approximated 

to the center-of-mass velocity at this order, Ve ■ The physical meaning of this equation is that electrons gain 

(or lose) their momentum through scatterings due to the relative velocity to photons, and the anisotropic stress of 
photons. The momentum transfer from the photons ensures the velocity difference between electrons and protons, and 
thus eventually generates magnetic fields. We found that the contribution from the curvature perturbation is always 
much smaller than that from the density contrast of photons. Eurthermore the tensor perturbation, i.e., primordial 
gravitational waves, is subdominant comparing with the scalar perturbation in the current observations (toI - It^ . 
Therefore, we have omitted the curvature perturbation and the tensor perturbation in Eq. (1451) when considering the 
evolution of magnetic fields. Equation (HSl) shows that the magnetic field cannot be generated in the first order. The 
right-hand side of Eq. (1451) contains two types of source terms, i.e., a purely second-order term and those that consist 
of the products of first order quantities. 

The first term in Eq. (H5|) is exactly the same as that discussed in [^. They have estimated the amplitude 
of magnetic fields from these terms by considering typical values at recombination. Here, we solve the equation 
numerically and obtain a robust prediction of the amplitude of magnetic fields in the standard ACDM cosmology. 


C. Scalar, vector, and tensor decomposition 

We devote this subsection to rewriting the evolution equation of magnetic fields in the context of the scalar, vector, 
and tensor decomposition approach. In the standard cosmological perturbation theory, the perturbations can be 
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decomposed into three modes, i.e., scalar, vector, and tensor modes [7^ as 

= wo(fe)of (fc)+ ^ ujxik)0l^\k) , 


A=±l 


Xzj{k) = Xiso{k)5ij + xo{k)0^f (k) + ^ xxik)0^,j\k) + ^ Xa{k)0^f{k) 


(46) 

(47) 


X=±l 


<t=±2 


where and are the orthogonal bases for Fourier modes with wave number k = kk and m = 0, ±1, and ±2 
represent the scalar, vector, and tensor modes, respectively. Note that, in our notations, the normalization is different 
from that adopted in Ref. Q. These inverses are given as 


wo(fc) = -Ol°\k)u;^{k) , u;x{k) = -O) ''''(fc)wi(fc) , 


(-A), 


Xo(fe) = ^0l°\k)xijik) , xxik) = -20b ''^(fe)xy(fc) , x^(fc) == |Ob ’">{k)xij{k) . 


(-A), 




(48) 

(49) 
(0) 


In this notation, magnetic fields can be decomposed into the scalar and vector modes, i.e., Bi = SqO) + 

Sa=±i BxOl^\ When we pull out the scalar mode from Eq. (1451) . we find that the right hand side of Eq. (l45l) 
vanishes, namely, 

dBo 


dt 


= 0 . 


(50) 


This is because magnetic fields consist of rotation of the vector potential, or in other words, magnetic fields do not 
have the scalar component. In contrast, the vector mode for Eq. (1451) is given by 


dBx{k) 

dt 


4CTTP7°^a 

3e 


(Afc) 




(Pki 

(27r)^ 


^bo(*l)4'Hfc2)t/^ri:A(fcl) 


(fcl )n7 0 (^2)3^1^;A (fcl, ^2 ) 


(51) 


where we define the function y^^^^(ki, k 2 ) as 




£i 4 ^ 
0 0 0 


h £2 e 

mi m 2 —m 


47r 


2^1 + 


47r 


. (52) 


where Y^^rnik) is the spherical harmonics. In the above equation, the multipoles should satisfy the condition that 
■^1 + ^2 + ^ = even because of a property of the Wigner-3j symbol. Furthermore, the triangle condition, k 2 = k — ki, 
is ensured by the delta function. Throughout this paper, we keep fc 2 in equations for simplicity of presentation even 
after being integrating out. This function obeys the following relations, 


3^,%^=(fei,fe2) = 3^,X(fe2,fci) , yl:^{kiM) = ^:^Yl^{ki)6,,,, , 3^;;:=(fci,fe2) = 


47r 


;0.^2/ 


47r 

2 £+ 1 


Ylra{k2)5i,i, . 

(53) 

To obtain the magnetic fields from Eq. m, we need to know the second-order relative velocity between baryons and 
photons for the vector mode. In the next section, we show the formulation of the second-order Einstein-Boltzmann 
system. 


III. SECOND-ORDER EINSTEIN-BOLTZMANN SYSTEM 

In this section, we derive the perturbed equations of the Einstein-Boltzmann system up to the second order. 
Throughout this paper, we work in the Poisson gauge of which the line element is 

ds^ = a^{rf) \—e^'^drj^ -\- 2uJidr]dx’' + -I- Xij) dx'^dx^] . (54) 

Under the Poisson gauge, the gauge conditions = x^^,j = 0 ^^nd the traceless condition x% = 0 are imposed on uji 
and Xij- We expand the metric perturbations as 41 = -|- -|- LUi = and Xij = 

where we neglect the first-order vector and tensor perturbations. The first-order vector perturbation is neglected in 
the standard cosmology because the vector perturbation has only a decaying mode in the first-order perturbation 
theory. The first-order tensor perturbation is observationally shown to be subdominant compared with the scalar one 
at first order 0. 
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A. Boltzmann equation 


From Eq. (I51L full computation of magnetic fields needs to evaluate the relative velocity between baryons and 
photons in the second-order perturbation theory. The velocity perturbations of the two fluids are described by the 
Boltzmann equation. 

The Boltzmann equation with binary collisions can be written as 

^(x^P'^) = C[/] , (55) 


where A is the affine parameter and C [f] is the collision term, i.e., in the case of the photon distribution function, 
the Thomson interaction between photons and electrons. Note that protons and electrons interact with each other 
through the Coulomb interaction. We can treat proton and electron fluids as a single component since these particles 
couple strongly by the Coulomb interaction. For the Boltzmann equation for dark matter or neutrinos, the collision 
term must vanish. 

To calculate the perturbed Boltzmann equation, it is useful to change the coordinate system from the Poisson gauge 
to the local inertial frame [d^. Since we consider the cosmological perturbations up to the second 

order, the distribution function is expanded as 

f{r],x,p,n) = + f^^'>{r),x,p,n) + ^f^^\r),x,p,n) , (56) 

where p and fi are the amplitude and the direction of the photon’s momentum, respectively. The zeroth-order 
distribution function, f^^\ri,p), is fixed to the Planck distribution. It is useful to define the brightness function which 
is given by 


( 77 , X, h) 


J dp p^f<^^’^\7j,x,p, n) 
J dp p^f(°'){r],p) 


where the denominator of the right-hand side is proportional to the mean energy density of photons. 
The angle dependence of the brightness function is expanded by the spherical harmonics as 


(57) 




I m— — 


(58) 


The coefficients are related to the density perturbation, velocity, and anisotropic stress for photons as Aq^q = 
aJ^q = and A^^q = respectively, as is shown Eqs. (I?T]) - (I551) . The Boltzmann equation of photons in 

terms of A) at first- and second-order is written as 


A 


( 1 . 2 ) 

l,m 


+ k 


C£+l,m A (1,2) 

2 £ + 3 ^+ 1 ."* 


A (1.2) _ q(1.2) 

2£ _ ^£,m ’ 


(59) 


where ci^m = 'JP — w?. A dot represents a derivative with respect to the conformal time rj and here we have 
translated from real space to Eourier space. The source term ^ can be expressed as 

Sfl{k, r,) = r,) + gfl{k, Ty) , (60) 

Here, C) ^ is the collision term that is proportional to fc, where fc is the differential optical depth which is defined by 

the number density of the electron Ug, scale factor a, and the Thomson scattering cross-section ctt as fc = —angaT, 
( 2 ) ( 2 ) 
and Gi ^ denotes the gravitational effects, i.e., the leasing and the redshift terms. The collision term C) ^ is related 

to Eq. (1^ as 




2£+l 

An 


Ylmin) 


f dp p^C^^Yr], k,p, fi) 
J dp p3f(°){r],p) 


(61) 


where C*-^^(7y, k,p, h) is the Fourier transformation of Eq. (I26|) . while the gravitational effect G^Y coming from the 
left-hand side of the Boltzmann equation (1551) with the same procedure as obtaining k,p, h). In this paper, we 
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call C(^m and Ge,m the scattering term and the gravitational term, respectively. The explicit forms of and 
are given in our previous study [tHI ]. 

The source terms of the first-order Boltzmann equation vanish when m ^ 0, because we consider only the scalar 
mode in the first-order perturbations. However, for the second-order perturbations, not only the scalar mode (m = 0), 
but also the vector (m = A) and tensor (m = a) modes arise due to nonlinear couplings, where A = ±1 and cr = ±2, 
respectively. 

In the case of massless neutrinos, one can set fc = 0 in the above equations because massless neutrinos interact 
with the other fluids only through gravity. We do not write down the hierarchical equation of neutrinos here since it 
is trivial. The distribution function of neutrinos is also expanded by the spherical harmonics. 


B. Tight coupling solutions of the vector mode 


To solve the second-order equations derived in Sec. IIII Al numerically, we should set up the initial condition of 
each perturbation variable. Thus we first solve the equations analytically with /cry <C 1 and using the tight coupling 
approximation, and hnd the initial condition at sufficiently early time for our numerical calculation. 

Deep in the radiation dominated era, photon and baryon fluids are tightly coupled because the opacity Ty, is large 
1^. Although the photon and baryon fluids would behave as a single fluid, there is a small difference 
in motion between photon and baryon fluids. For this reason, we can expand the perturbation variables using the 
tight-coupling parameter which is given by 


, 10-2 ( _ - _^ 

I 1 A/Tvv^-i \ n no / 


-1 


(62) 


where Db is the baryon density normalized by the critical density, and h is the normalized Hubble constant. In what 
follows, we derive the tight-coupling solution up to the first order to set the initial condition of photon and baryon 
fluids at second-order in cosmological perturbations and to calculate the evolution of perturbations in a numerically 
stable manner. 

We expand the cosmological perturbation variables using the tight-coupling parameter e up to the first order as 

^(CPT=l,2) ^ ^(CPT=l,2,TCA=0) ^(CPT=1,2,TCA=/) _ /ggN 


where the Arabic number and the Roman number represent orders in the cosmological perturbation and the tight 
coupling expansion, respectively. In the rest of this paper, we focus on the vector mode m = ±1. 

First, the solutions at zeroth order in the tight-coupling expansion, namely, in the tight-coupling limit, are given as 




0 , 


20 


(Pki 


■^7 0 


y2x{ki,k2) , 


(64) 

(65) 

( 66 ) 


where k 2 = k — ki in Eq. (I65|) . In the tight-coupling limit, the relative velocity between photons and baryons vanishes 
as well as in the first-order cosmological perturbation theory. However, the anisotropic stress of photons is present 
due to the quadratic of the photon velocity. 

Second, the solutions at first-order in the tight-coupling expansion are given as 
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1 + r (2,/) 

- ov. 


R 


7b A V 


(fc) = ^ - (^) + 


+ 


(Pki 
d?ki 
(Pki 
f <Pki 




7 b 0 

i^^b^o(^l)^^y^o(^2) 


24'o(fci)Kb'o(fc2) 


7 

,( 1 ) 


-^7(^1) 


i? 


+ 


/ (27r)3 
(fiki f ki 




(27r)3 \ 
(Pki 


-i<5W(fci)(vI/(i) + <i>(i))(fc2) - 2vI/(i)(fci)5W(fc2) 


+ 

+ 


(27r)3 

(Pki 
(27r)3 

d^fci / /C2 
(27r)3 \ T, 


1 

-V. 


-.W(fci)(<5W+4vI/W) (fc2) 


/ 47r 


-3^7^o(^i)'^7^^(^2) 


- 3 ^^ 071 ) 4 ^ 072 ) 


/ 47r 

y777i) 


-77 7i) 


1 


(Pki 

WF 

f (fiki 


-^n4o7i)4472) 


747i72) 


/ (2^)3 
d3fc]^ / ki 


-^nV7fci)44(fc2) 


(2^)3 


7f(i) 
2 

3 


yliikiM) 


-7-77(7)4^47) 


3^i7(fci,fe2) 


- 77 ( 7 ) 


(67) 


4-A477fc) = -f 1^43 A (2) 

10 2.A 1 1 \f 


I 


d?ki 

dPki 


Al 7 (fe) 

-9nW(fci)(<5W+vI/W)(fe) 


/ 47r 

'y77(7) 


(27r)3 \^rc. 




3 


/ 47r 

'y>^27(7) 


184^o(7)4^o( 7) + 84^o(7)<54b o(7)j 3^24(777 



)-8$(i)) (fci)47(fc2)] 3^24(7,7) 


a(2.^) 

^3,A 


7 ) = - 




cP'ki ' 

Wf 


< 47 ) 44 ( 7 ) 13 ^ 34 ( 772 ) 


Af74fc) = 0 


( 68 ) 

(69) 

(70) 


where R = 3p\^ 7(47 ). To derive the above solutions, the time derivative of the first-order curvature perturbation 
$ is ignored due to the conservation of the curvature perturbation on large scales. In the right-hand side of Eqs. (EH), 
(IB51) . and (ESI), for simplicity, we omit the superscript for the order of the tight-coupling parameter. Note that the 
octopole does not vanish and higher multipoles than i = 3 are equal to zero at this order [75l| . 
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C. Einstein eqnation of the vector mode 


In this subsection, we derive the second-order Einstein equation and hereafter we focus on the vector mode only. The 
evolution equation for the metric perturbation of the vector mode, i.e., uj\ can be derived from the (*, j) component 
of the Einstein equation as 

a^G^j = -f 

-|-(diagonal part) S^j , (71) 


and 


T^j = PiUlj + (diagonal part) , (72) 

T"n,j = + (diagonal part) 5^j , (73) 

where T^j and T^j represent the energy-momentum tensors of massless (relativistic) particles such as photons and 
neutrinos, and massive (nonrelativistic) particles such as baryons and dark matter, respectively. The equation of the 
vector mode is derived by acting the projection operator, —0\^ which can pull the vector mode, as 




2 1 


5^3 fc 


(87rGa2p^°)Ag(fc) -f SirGa'^p^°'>J\f^‘^l{k)'^ 


(Pki 


4fci 


(27r)3 

cPki 4 fcf r 
(27r)3 .^3 k I 


■ / /[ 


1 I zL-tt 

$W(fci)$W(A:2) -f T(i)(fci)«'W(fc2) \ 


+ X! 87rGaV(°^ 


s=b,dm 


(Pki 

Wf 




(1)/ 


YKoi^i)\jY^i:xik2) . (74) 


Note that in the Poisson gauge, the vector metric perturbation is only included in uji. When we consider the evolution 
equation up to the first order in the standard cosmology with perfect fluids, the right-hand side of Eq. (HI) becomes 
zero. As a result, the vector mode has only a decaying solution, which is neglected in the linear theory. 


IV. COSMOLOGICAL MAGNETIC EIELDS 

In this section, we show the evolutions and spectra of magnetic fields driven by the Harrison mechanism. In the 
previous studies 0 m m, [zi, generated magnetic fields are partially estimated by numerical or aiialytical 
ways, and there is a small discrepancy between in Refs. and [i^. In this paper, we build on Ref. and 
expand the work by including all contributions numerically. The source terms of the magnetic fields consist of three 
contributions, i.e., and which hereafter we call “the slip term”, “the anisotropic stress term” 

and “the second-order slip term”, respectively. 

In this paper, we focus on three issues on the generation of magnetic fields at recombination. First, we consider 
how large is the contribution of the second-order slip term on magnetic fields compared with the contributions of the 
slip and the anisotropic stress terms. Second, to evaluate the total spectrum of magnetic fields we need to include 

the cross-correlation terms between the sources, namely, Pb ^ (^{B 2 nd Slip + -Bsiip + SAnis)^^ The cross terms can 
be negative and it has the possibility to cancel the generated magnetic fields from each of the source terms. Third, 
we try to find the cause of the small discrepancy between Refs. [63 and [gj. In the following subsections, we show 
the evolutions and spectra of magnetic fields and the answers of the above three considerable questions. 


A. Configuration in Fourier space 

Before moving to the results, we mention the configuration of Fourier space in carrying out the convolution since all 
second-order quantities are written by the products of the first-order scalar perturbations. To calculate the second- 
order power spectrum, we decompose the second-order variable into the transfer function and the primordial amplitude 
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M 


as 


A(2)(?7,fc) = y J ^^6{k-ki-k2)A^^\r],k,ki,k2)^{ki)^{k2) , (75) 

f2l 

where and $ are the second-order transfer function and the primordial amplitude, respectively. The ensemble 
average of the variance of the primordial amplitude can be expressed as ($*(fci)$(fc 2 )) = {2TT)^P^{ki)S{ki — ^ 2 ), 
where is the primordial power spectrum determined from cosmological observations such as CMB and large 

scales structure. In this paper, we use the power-law spectrum as 


4 

= 5 A 



(76) 


where the parameters. As and rig, are the amplitude and the spectral index of primordial perturbations, respectively. 
We set As = 2.4 x 10“® from the WMAP nine-year results [s^, and for simplicity, we consider a scale-invariant 
spectrum, namely, rig = 1.0. 

In this decomposition, the second-order power spectrum can be written as 


P^(2){r],k) = 2 


d?ki 


A^^\r],k,ki,k2) 


P^ik\)P^{k2) , 


(77) 


where we can use the symmetry under the exchange of ki and ^2 without loss of generality to derive the above 
relation. Note that k 2 = k — ki should be satisfied implicitly, namely, A^\ri, k, ki, ^ 2 ) = A^\rj, fc, fei, k — ki) and 
T$(^ 2 ) = P<s>{\k — fell), as is mentioned before. 

We need to solve the Einstein-Boltzmann system in {k, fci, ^ 2 ) space. Note that the transfer function is transformed 
under the rotation of (j) as 


A^")(77,fc,fci,fc2) ^ A^")(77,fc,fci,fc2)e™‘^ . (78) 

In practice, we take (j) = 0 = 0^ (f)i = 0^ and (j )2 = for k, fei, and k 2 , respectively. In other words, the transfer 
function under the exchange of fci and ^2 transforms as 

A^")(7?,fc,fci,fc2) = (-l)™A^2)(7?,fe,fc2,fci) . (79) 

It is very interesting that for the case of the scalar and tensor modes, the dominant contribution comes from near 
the fc ~ fci ~ ^ 2 . Because the transfer functions of the scalar {m = 0) and tensor (m = ±2) modes are symmetric 
under the exchange of ki and k 2 - On the other hand, the vector mode does not have dominant contribution near the 
fc ^ fci ~ ^2 since the transfer function of the vector mode is antisymmetric under the exchange of fci and k 2 - 

In the vector mode, we can naively consider three configurations of the triangle k = ki + k 2 that contribute to the 
power spectrum on superhorizon scales, namely, the one where both ki and ^2 are on superhorizon scales, where ki 
(^ 2 ) is at superhorizon scales while k 2 (fci) is at subhorizon scales, and where both ki and ^2 are at subhorizon scales. 
When both ki and k 2 are at subhorizon scales, the square of the primordial power spectrum in Eq. (Eli) does not 
contribute on the power spectrum because of the antisymmetric nature of the vector transfer function. In the vector 
mode, therefore, we need a careful treatment of how to sample wave numbers in ik, fci, ^ 2 ) space. This difficulty does 
not arise in the second-order scalar and tensor modes calculations because the transfer functions of the scalar and 
tensor modes have a symmetry under the exchange of fci and ^ 2 . 


B. Evolutions of magnetic fields 

In Fig.[T] we show the evolutions of magnetic fields induced by the slip, the anisotropic stress, and the second-order 
slip terms. First, we focus on the evolutions of magnetic fields before the horizon crossing. The time evolutions of 
magnetic fields can be approximated by a power law, and the powers of the slip and the anisotropic stress terms are 
T^siip oc ® and BAnis oc respectively. These results correspond to the previous study [i^. Furthermore, the 
second-order slip term is proportional to rf'^ on superhorizon scales. We can explain the coincidence of the powers 
between the anisotropic stress and the second-order slip terms as follows. In Eq. dSZl), the dominant terms in early 
times can be estimated by using the superhorizon solutions at linear order. At first glance, the term x 
in Eq. (1571) seems to give a dominant contribution to the second-order slip term. However, because both and 
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FIG. 1: Evolutions of generated magnetic fields sourced by the slip term {top left), the anisotropic stress term {top right), and 
the second-order slip term {bottom). We show evolutions of generated magnetic fields for wavenumbers k = 10~^ hMpc“^, 
10“^ /iMpc“^, and 10° /iMpc“^ as indicated in the above panels. We can see that magnetic fields at smaller scales generated 
earlier and their amplitudes are larger. 


are scale invariant in the Poisson gauge and according to the formula of spherical harmonics, y + 

= kSm,o, this scale invariant term vanishes in the vector mode. As a result, we find that the most 

dominant term in Eq. (EZl) is that proportional to which is the same form as the anisotropic stress term. 

Therefore, the powers of the time evolutions of the anisotropic stress and second-order slip terms coincide. 

Next, we discuss the evolutions of magnetic fields after the horizon crossing. We can see that the magnetic fields 
from the slip and the anisotropic stress terms start to decay adiabatically as oc a“^, after the horizon crossing since 
their sources also diminish after the horizon crossing. On the other hand, magnetic fields induced by the second-order 
slip term do not decay adiabatically even after the horizon crossing. This arises from the fact that the second-order 
slip term continues to grow even after the source terms from the first-order perturbations become negligible, 

until the corresponding scale reaches the Silk damping scale. Therefore, the purely second-order perturbations can 
contribute to the magnetic field generation even if the product of the first-order perturbations is absent. When we 
neglect the product of the first-order perturbations, the evolution equations for magnetic fields are corresponding to 
the case of the first-order magnetic fields generation [s^ . The relative velocity between photons and baryons, ^ 
contributes to the generation of magnetic fields after the horizon crossing. This additional enhancement can be seen 
in the bottom of Fig.[T]at k = 10° /iMpc”^. However, magnetic fields induced by the additional enhancement undergo 
nontrivial cancellation after the Silk damping epoch and magnetic fields decay faster than the adiabatic decay that it 
is proportional to [s^. The final amplitude of magnetic fields is consequently determined by the initial amplitude 
around horizon crossing. 
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k [/zMpc'^] 


FIG. 2: Spectra of generated magnetic fields sourced by the slip term {top left), the anisotropic stress term {top right), and 
the second-order slip term {bottom). We show spectra of generated magnetic fields for redshifts z = 1.1 x 10®, 1.1 x 10^, and 
1.1 X 10® as indicated in the above panels. 


C. Spectra of magnetic fields 


Next, we show the spectra of magnetic fields induced by the slip, the anisotropic stress, and the second-order slip 
terms in Fig. O From Fig. [21 resultant magnetic fields are dominated by the anisotropic stress and second-order 
slip terms on smaller scales, i.e., k > 1.0 /iMpc”^, at 1 + z = 1100. Conversely, on these scales, the slip term is a 
subdominant source for magnetic fields. 

On superhorizon scales, we can see that the spectra of magnetic fields are proportional to /c® ®, which also corresponds 
to the results about the slip and the anisotropic stress terms in Ref. . This power is also consistent with the power 
spectrum for causal magnetic fields [8l| . From Ref. [6^ , the magnetic power on superhorizon scales can be estimated 
as below. For example, we focus on magnetic fields induced by the slip term. We can integrate the evolution equation 
for the second-order magnetic fields (l?T]) and take the ensemble average as 


27r^ 


Psik) 


oc k° 


Slip 


d?ki 


(l - pil) P<5(fcl)P$(fc2) 


S\k2,ki) -^S{ki,k2)S{k2,ki) 

k2 


where fii = cos0i and S{ki,k 2 ) is defined as 


(80) 


S{k^,k2) 


J dr]a^{'n)p^°\r])d!:^^^ki,r])dv^^2oi>^2,r]) , 


(81) 


where for simplicity, we omit the time dependence of S{ki, ^ 2 )- To proceed the estimation of the magnetic power, we 
take the limit k/ki —>■ 0. This approximation can include the contributions from subhorizon scales. In this limit, 


1 - + O [{k/ki)"^) 

ki 


k2 = fci 


(82) 
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Furthermore, we can approximate the integrated source term as S'(fci, ^ 2 ) ~ S{k 2 , fci) ~ T{ki) since S{ki, ^ 2 ) can be 
treated as k independent in the above limit. Then by using the fact that oc Eq. (IM)) can be rewritten as 


27r2 


PB{k) 


Slip 


OC 


k^ J kldki J dfii (1 - fil) k^-‘-*k^-‘-* (^1 - T2(fci) 

k^ . 


(83) 

(84) 


This nonlinear power law can be seen in Fig. [2] Note that if we use the superhorizon solution only, namely, oc (kij)^ 
and Q (X k^r/^, the magnetic power is returned as oc fc®. This superhorizon power does not match for our numerical 
results. 

By using our numerical code, we trace a possible cause of the discrepancies between the results in previous studies. 
As we noted in Sec. IIV A[ the transfer function of the vector mode is antisymmetric under the exchange of fci and 
^ 2 , and therefore, the isosceles conhguration such that ki = ^2 in Fourier space does not contribute in the calculation 
of the power spectrum. However, to achieve the result correctly, contributions from the configurations of ki < ^2 
and ^2 ^ ki should be included. When these contributions are not included in the numerical calculation, the power 
spectrum of magnetic fields on super horizon scales shows oc which corresponds to the result obtained in Ref. [b^. 

On subhorizon scales, we can find that the spectra of magnetic fields induced by the slip, the anisotropic stress, 
and the second-order slip terms are proportional to and respectively. The spectra induced by the 

slip and the anisotropic stress terms are consistent with Ref. [63|. A noticeable feature in the spectrum induced by 
the second-order slip term is the additional amplihcation at fc « 5.0 x 10“^ /iMpc”^ at z = 1100. As discussed 
in Sec. IIVBI this additional amplification is due to the second-order relative velocity between photons and baryons 
after the horizon crossing. However, this amplihcation is a temporary effect and the amplihed magnetic helds by 
this effect had been erased by the epoch of Silk damping [s^. Therefore, the amplihcation cannot be seen for scales 
k > 1.0 /iMpc"^ 

The above discussions and results in this subsection are valid only for the auto-power spectra of the magnetic helds 
from the slip, the anisotropic stress, and the second-order slip terms. In the following subsection, we focus on the 
total power spectrum induced by all the contributions including the cross spectra. 


D. The second-order magnetic fields 


We depict the total power spectrum at recombination in Fig. [3l It is clear that the second-order slip term gives 
a dominant contribution to the total magnetic fields. The amplitude of magnetic fields from the second-order slip 
term is 10 times larger than without the second-order slip term. However, we find that the magnetic fields from the 
second-order slip term are canceled out by the magnetic fields from the anisotropic stress term on small scales. 

By using the tight coupling solution given by Eq. (Ell), this cancellation is easily understood analytically. As we 
mentioned in Sec. IIVBI the dominant term in Eq. (1671) is the anisotropic stress term given by 


5v 


(2) _ 
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\ + R 


d^ki 
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-^n^^o(b)w[^o(fc2) 


yiiiki.h) 


(85) 


while any other terms related to the anisotropic stress are subdominant with the baryon-photon ratio R oc 

as a suppression factor. By substituting this expression into Eq. dm), the evolution of magnetic fields from the second- 

order slip and the anisotropic stress terms is given as 


dBx [5 1 

dt 4 1 -|- i? 
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'^b 0 


ik,)n%{k2)yl:lik,k) , 


( 86 ) 


where the first and the last terms in the parentheses are coming from the second-order slip term and the anisotropic 
stress term, respectively. We can see that the two terms in the parentheses in Eq. (1551) are canceled in the radiation 
dominated era, where R is negligibly small. However, in the matter dominated era, the baryon-photon ratio has large 
value and this cancellation does not occur. 

The dominant contribution from the second-order slip term is canceled out by the contribution from the anisotropic 
stress term in the radiation dominated era. Conversely, there still remain some contributions from the second-order 
slip term as shown in Eig.j^and discussed below using the tight-coupling solution. The subleading contribution from 
the second-order slip term can be written as 
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(87) 
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FIG. 3: Magnetic spectra generated from the slip term, the anisotropic stress term, the second-order slip term, and all terms 
included the cross terms at recombination (1 -h z ~ 1100). 


Then, the evolution equation of magnetic fields induced by the slip and the second-order slip terms can be rewritten 
as 


Zx(k{) , ( 88 ) 

where the first and the last terms in the parentheses are coming from the second-order slip term and the slip term, 
respectively. From Eq. (l88)) . we find that the total amplitude of the spectrum of magnetic fields is twice as large as 
the case only with the slip term. This tendency can be seen in Fig. [31 

Next, let us discuss the evolution of magnetic fields through the epoch of recombination. As the process of re¬ 
combination proceeds, electrons form neutral hydrogen atoms with protons and the number of free electrons rapidly 
decreases. Accordingly the effect of the Compton scattering on generation of magnetic helds becomes negligible and 
magnetic fields are no longer generated through the Harrison mechanism. We show the power spectrum of the second- 
order magnetic fields after recombination at z ~ 500 in Fig. [T] In Fig. [31 one can find new features in the spectrum 
different from one at recombination on intermediate scales such as 10“^ ft,Mpc“^ ^ k < 1.0 ZiMpc”^. These features 
are nearly consistent with Ref. [^ . For instance, magnetic fields induced by the slip term decrease on intermediate 
scales. On the other hand, magnetic fields induced by the second-order slip term are enhanced at cosmological re¬ 
combination since the relative velocity is also enhanced at that era. We furthermore extend the magnetic spectrum 
to much smaller scales. On scales where 1.0 ft,Mpc“^ > fc, magnetic fields induced by the slip, the anisotropic stress, 
and the second-order slip terms have the same structures as the magnetic spectrum at recombination shown in Fig. (3] 
As a result, the spectrum of second-order magnetic fields has a slightly blue tilt on small scales. 


dBx 

dt 


1 


1 + R 




E. Nonhelical magnetic fields 

Before closing this section, we investigate the possibility whether helical magnetic fields are generated in the Harrison 
mechanism or not since the helicity may play important roles in the cosmological observations In fact, in 

Ref. [13, the authors found the evidence of existence of helical magnetic fields on a few Mpc scales. It is believed that 
helical magnetic fields can only be generated through the process of parity violation. Because helicity is conserved 
in the standard magnetohydrodynamics, it is a good indicator to probe the generation mechanism of magnetic fields. 
We will show below that the Harrison mechanism does not induce helical magnetic fields since this mechanism relies 
on the standard Compton scattering which does not break the parity symmetry. 
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FIG. 4: Same as Fig. [3] but for after recombination as (1 + ^ ~ 500). 


At first, under the existence of helical magnetic fields, the correlation of magnetic helds can be written as 


{B,{k)B,{k')) 


{2'Kf5{k - k') 




4“ 


Pnik) 

2 


(89) 


where Psik) and Pnik) are the spectra of nonhelical and helical magnetic fields, respectively. The power spectrum of 
helical magnetic fields can be pulled by the subtraction as {Bi{k)Bj{k')) — {Bj{k)Bi{k')). In the Harrison mechanism 
that is given by Eq. dill), generated magnetic fields can be symbolically expressed as 


B^{k) oc Ciabk^" 


(Pki 


k\f{k^M)X^^\k^)Y'^^\k2) , 


(90) 


where /(fci,fc 2 ) is an arbitrary real function of fci and ^ 2 , and X^^'>{ki) and Y^^\k 2 ) are the time integrals of first- 
order scalar perturbations. The scalar perturbations can be decomposed into the primordial perturbation ^^P{k) 
and the transfer function AT(fci) or TT(fc 2 ) as X^^'>{ki)Y^^'>{k 2 ) = (fc 2 )AT(fci)FT(fc 2 )- Note that we 

use the fact that the purely second-order variables, e.g., Sv^ti, are composed the product of the first-order scalar 
perturbations. 

Finally, we evaluate the helical part of the power spectrum as 


{Bi{k)Bj{k')) - {Bj{k)Bi{k')) oc {2Tr)^6{k - k') {eiateja'b' - ijab^ia'b') k°-k° 


(Pki 

(^ 


/(fci, fc2)^T(fci)Tr(A:2)H$(fci)P$(fc2) 


X 

= 0 


fcJfcf/(fei,fe2)XT(fci)yT(fc2) + i (k\ki+klk^;) /(fc2,fcl)^T(fc2)bT(fci) 


(91) 


where we symmetrize about 5 o 6' in the square bracket by using the nature of the symmetry under the exchange of 
ki and k 2 - In conclusion, the Harrison mechanism cannot induce helical magnetic fields. This result is coming from 
the fact that general relativity and the standard Maxwell theory do not violate the parity symmetry. Therefore, the 
observed helical magnetic fields call for other mechanisms to explain. 
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V. CONCLUSION AND SUMMARY 


In this paper, we reinvestigate the spectrum of magnetic fields induced by cosmological perturbations through 
the Harrison mechanism. If we consider the cosmological perturbation theory up to the first order, the Harrison 
mechanism does not work since the vector mode, which is needed for this mechanism, has only a decaying solution. 
However, when we expand the cosmological perturbations up to the second order, the regular solution of the vector 
mode is excited by the first-order scalar mode. The Harrison mechanism works in the higher-order cosmological 
perturbation theory. 

In previous studies, the spectrum of magnetic fields induced by this mechanism has been estimated. In Ref. the 
authors show the spectrum of second-order magnetic fields induce d by the product of the first-order perturbations, 
namely, the slip and anisotropic stress terms. Subsequently, in Ref. [6^, the purely second-order slip term is included. 
By comparing these works, however, it is found that there are some discrepancies in the product of first-order 
perturbations. For example, the power law tails of the spectrum induced by the slip and anisotropic stress terms on 
large scales have different fc-dependences in Refs, and [gj. Furthermore, the scale-dependences of the spectrum 
are slightly different from each other. We find that the discrepancy can be explained by the lack of sampling in the 
Fourier modes at fci ~ ^2 of the first-order scaler perturbations in Ref. [g^, and our results agree with the ones of 
Ref. [H. 

Let us summarize features of the magnetic fields induced by the second-order magnetic fields at cosmological 
recombination as follows. 

• The scale dependence of magnetic fields on large scales is oc which is consistent with the result in Ref. [g^. 
Note that magnetic fields generated by causal processes have the same power [8l| . 

• On small scales, the spectra of magnetic fields induced by the slip, the anisotropic stress, and the second-order 

slip terms have the power of and k^-^, respectively. In particular, the spectra of magnetic fields 

induced by the slip and anisotropic stress terms are consistent with the result in Ref. [gj- 

• The cancellation occurs between the anisotropic stress term and the second-order slip term on small scales in 
the tight coupling regime in the radiation dominated era, and the power of magnetic field spectrum becomes 
oc This result indicates that the spectrum of magnetic fields cannot have the large amplitude as argued in 

Ref. [13. 

• The spectrum of magnetic fields at cosmological recombination has a bump at fc ~ 5.0 x 10“^ /iMpc”^ owing to 
extra amplihcation after the horizon crossing, where the amplitude of magnetic fields is Spec ~ 5.0 x 10“^"^ Gauss. 
However, after all, this amplification vanishes by nontrivial flipping of the relative velocity between photons and 
baryons discussed in Ref. [33 ■ This cancellation creates characteristic diffusion scales in the magnetic fields 
spectrum around k r; 10° /iMpc”^. 

• The Harrison mechanism does not work efficiency below the scale of Silk damping at the electron-positron pair 
creation epoch as fccut ~ 10° liMpc”^, as discussed in Ref. [^. Even if we extrapolate our numerical result 
toward smaller scales by using the power of oc ° up to the cutoff scale, the amplitude of the magnetic fields at 
that scale cannot be larger than 10~^° Gauss, assuming that the linear density perturbations are scale invariant. 


Finally, we discuss implications of the cosmological seed fields. The derived amplitude of magnetic fields at recom¬ 
bination has a peak about 5.0 x 10“^^ Gauss, which is sufficient for a candidate of the seed of galactic magnetic fields 
However, this amplitude seems to be somewhat small to explain the intergalactic magnetic fields [13, [ll|. Note 
that, the above amplitude is derived assuming that the primordial perturbations are scale-invariant, while primordial 
perturbations with a blue tilt lead a larger amplitude of the magnetic fields on smaller scales. 

The magnetic fields induced by the second-order perturbation must be inevitably generated in the standard cos¬ 
mology, and it is possibly that the magnetic fields act as seed fields for the turbulent dynamo during the structure 
formation of the universe. 
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